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Abstract 

Phenomenologically viable and interesting regions of parameter space in the minimal 
super-gravity (mSUGRA) model with small mo and small mi/2 consistent with the 
WMAP data on dark matter relic density and the bound on the mass of the lightest 
Higgs scalar > 114 GeV from LEP2 open up if the rather adhoc assumption Ao=0, 
where is the common trilinear soft breaking parameter, employed in most of the 
existing analyses is relaxed. Since this region corresponds to relatively light squarks 
and gluinos which are likely to be probed extensively in the very early stages of the LHC 
experiments, the consequences of moderate or large negative values of are examined 
in detail. We find that in this region several processes including lightest supersymmetric 
particle (LSP) pair annihilation, LSP - lighter tau slepton (fi) coannihilation and 
LSP - lighter top squark {ti) coannihilation contribute to the observed dark matter 
relic density. The possibility that a ti that can participate in coannihilation with the 
lightest neutralino to satisfy the WMAP bound on relic density and at the same time be 
observed at the current experiments at the Tevatron is wide open. At the LHC a large 
number of squark - gluino events lead to a very distinctive semi-inclusive signature 
t^+Xt- (anything without a tau lepton) with a characteristic size much larger than 
e^+y^e or /i^-|-X^ events. 
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1 Introduction 



Models with supersymmetry(SUSY)[l] are interesting for a variety of reasons. Theoretically 
the removal of quadratic divergence in the Higgs boson mass in the standard model (SM) 
by similar divergent loop diagrams involving SUSY particles (sparticles) is very attractive. 
These models have also predicted very interesting experimental signatures and are very 
promising candidates for beyond the standard model physics. A specially attractive feature 
of the minimal supersymmetric standard model (MSSM) with R-parity conservation is the 
presence of the stable, weakly interacting lightest neutralino {xi ) [2] which turns out to be 
a very good candidate for the observed cold dark matter (CDM) in the universe. [3, 4, 5]. 

Since superpartners are yet to be observed supersymmetry must be a broken symmetry 
and one requires a soft SUSY breaking mechanism that preserves gauge invariance and does 
not give rise to any quadratic divergence. The MSSM has a large number of such soft 
breaking parameters. Being guided by well motivated theoretical ideas as well as low energy 
phenomenology, models with specified SUSY breaking mechanisms drastically reduces the 
large number of such arbitrary parameters to only a few. The N = 1 supergravity models 
incorporate gravity mediated supersymmetry breaking and the models are attractive because 
of many features like gauge coupling unification, radiative breaking of electroweak symmetry 
[6] , controlling flavour changing neutral currents (FCNC) by specific and simple assumptions 
at the unification scale [1]. 

The simplest gravity mediated SUSY breaking model - the minimal supergravity 
(mSUGRA)[7] model has only five free parameters. These are the common gaugino and 
scalar mass parameters mi/2 and mo, the common tri-linear coupling parameter Aq, all given 
at the gauge coupling unification scale {Mq ~ 2 x 10^^ GeV), the ratio of Higgs vacuum 
expectation values at the electroweak scale namely tan j3 which is in fact largely independent 
of scale and the sign of fi, the higgsino mixing parameter. The magnitude of fi is obtained 
by the radiative electroweak symmetry breaking (REWSB) mechanism[6]. The low energy 
sparticle spectra and couplings at the electroweak scale are generated by renormalization 
group evolutions (RGE) of soft breaking masses and coupling parameters[8, 9]. 

Various SUSY models have been constrained by the data on cold dark matter relic 
density. But the recent resurgence of interest in this field is due to the very restrictive data 
from the Wilkinson Microwave Anisotropy Probe (WMAP) observation [10]. Combining the 
WMAP data with the results from the SDSS (Sloan Digital Sky Survey) one obtains the 



1 



conservative 3 a limits 



0.09 < VtcDMh^ < 0.13 (1) 

where VLcDuh^ is is the DM rehc density in units of the critical density, h = 0.71 ± 0.026 
is the Hubble constant in units of 100 Km s~^ Mpc~^. In supergravity type of models Xi 
becomes the LSP for most of the parameter space[4, 5] and one may consider Qcdm = ^x?' 
We should note here that the upper bound of Q^q in Eq.(l) is a strong limit but lower bound 
becomes weaker if we accept other candidates of dark matter. 

In the thermally generated dark matter scenario, at very high temperature of the 
early universe (T >> Tn.^o), Xi ^^s in thermal equilibrium with its annihilation prod- 
ucts. The annihilation products include fermion pairs (//), gauge boson pairs (iy"'"iy~ 
& ZZ), Higgs boson pairs [hh, HH, AA, hH, hA, HA, H'^H~) or gauge boson-Higgs boson 
pairs {Zh,ZH,ZA & VT^iy^) through s,t and u channel diagrams. Thereafter, at lower 
temperatures the annihilation rate falls below the expansion rate of the universe and Xi 
goes away from thermal equilibrium (freeze-out). The present value of the Q^oh^ can thus 
be computed by solving the Boltzmann equation for n^o, the number density of the LSP in 
a Friedmann-Robertson- Walker universe. Finding the neutralino relic density most impor- 
tantly involves computing the thermally averaged quantity < cTeffV >, where v is the relative 
velocity between two neutralinos annihilating each other and (Tg// is the neutralino annihi- 
lation cross section which involves all possible final states. In addition to the annihilation 
of a LSP pair, coannihilations of the LSP [11, 12, 13, 14, 15] may also be important. This 
happens if there are sparticles with masses approximately degenerate with the LSP mass. 

The annihilation cross section ae/f depends on the magnitude of the bino (5), the wino 
(W) and the Higgsino {Hf , H2) compositions of the lightest neutralino and proximity in 
mass of the LSP with any coannihilating sparticle. Generically the LSP composition is given 
by the following mixing: 

X? = NuB + N,2m + N,,H° + N,,Hl (2) 

Here the coefficients Nij are the elements of the matrix that diagonalizes the neutralino 
mass matrix. One typically quantifies the composition through the gaugino fraction of Xi 
which is defined as Fg = {Nul"^ + |A'^i2p. A xi would be called gaugino like if Fg is very 
close to 1(~ 0.9) , higgsino like if Fg ~ 0.1. Otherwise the LSP would be identified as a 
gaugino-higgsino mixed state. 
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Accordingly a typical MSSM parameter space (with gaugino mass universality) where 
the WMAP constraint is satisfied can be classified into several regions depending on the LSP 
annihilation/coannihilation mechanism. The list goes as follows. 

i) In the so called stau coannihilation region, a large degree of Xi — coannihilation[ll] 
reduces the thermal abundance sufficiently so as to satisfy the WMAP data. The WMAP 
allowed region in the mo— mi/2 plane in the mSUGRA model typically falls near the boundary 
of the forbidden parameter space where fi becomes the LSP. The coannihilation processes 
are of the type x?4 ^ ^"7. 44 ^ and 44 ^ ^'"^^ 77, 7^, ZZ, W+W- , hh. 
Here / is essentially the fi. 

ii) The focus pomt[lQ\ or the Hyperbolic hranch[ll] region in the mSUGRA model is 
characterized by a reduced value of A small causes the LSP to have a signif- 
icant Higgsino component or it can even be a pure Higgsino. Strong coannihilation of 
the LSP with lighter chargino occurs in this zone. Some of the dominant coannihi- 
lation processes in these region are[13, 14]: XiX^7X2Xi' ~^ Uidi,eih'i, AW~^ , ZW^ ,W^h; 
xt Xi 7 X1X2 UiUi, didi,W~^W~ . Having the smallest mass difference between coannihilat- 
ing sparticles the process XiXt indeed dominates among the above channels. As a result the 
thermal abundance is reduced appreciably so that it satisfies the WMAP data or coannihi- 
lations may even reduce it further (below the lower limit of Eq.l) thus causing Xi to be a 
sub dominant component of dark matter. 

iii) The funnel or the Higgs-pole region[18, 19] satisfies the WMAP data for large values 
of tan/3 extending to both large mo and large mi/2 regions. This is characterized by the 
direct-channel pole 2m^o ~ m^,m/^. 

iv) For a limited range of large negative values of Ao, the lighter stop can become very 
light so that it may coannihilate with the LSP. This is the stop coannihilation region[12] 
characterized by a very light m^^. 

v) In the bulk annihilation region or the bulk region[ll] where mo and mi/2 are reason- 
ably small in mSUGRA, the LSP turns out to be bino dominated and, consequently, couples 
favourably to right sleptons, which in fact are the lightest sfermions in this region of param- 
eter space. As a result an LSP pair annihilates into SM fermions via the exchange of light 
sfermions in the t-channel. This cross section depends on the mass of the LSP (jn~o), its 
coupling with the sfermions and the masses of the exchanged sfermions [18, 4, 5]. However, 
LEP2 has imposed strong bounds on sparticle masses[20], particularly on the slepton masses 
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in the present context of mSUGRA and typical studies with Aq = disfavors a part of the 
bulk annihilation zone. Additionally, a very severe restriction appears on the (mo — mi/2) 
plane of mSUGRA from the bound on lightest Higgs boson mass {rrih) [21] which practically 
rules out the entire annihilation region. For rrih >114 GeV one finds that within the frame 
work of mSUGRA the sleptons are significantly heavier than the direct LEP2 bounds and 
this leads to very small LSP annihilation cross section implying an unacceptably large relic 
density. Thus it has often been claimed in the recent literature that the mSUGRA param- 
eter space with low values of both mo and mi/2 is strongly disfavoured. This automatically 
eliminates the bulk annihilation region [22]. 

For small mo but relatively large mi/2 one obtains a parameter space consistent with the 
Higgs mass bound. In this region a large mi/2 pushes up the mass of the lightest neutralino 
while the lighter tau slepton (fi) mass is rather modest because of a small mo. Thus the LSP 
and fi are approximately mass degenerate and LSP-fi coannihilation [11] provides a viable 
mechanism of producing an acceptable relic density. However, the relatively large mi/2 tends 
to push up the masses of squarks and gluino and this reduces the size of the LHC signals 
significantly. 

The purpose of this paper is to emphasize that the above conclusions are artifacts of the 
adhoc choice = 0. Many of the current analyses invoke this choice without any compelling 
theoretical or experimental reason. On the other hand it is well known that moderate to 
large negative^ values of Aq lead to larger m^[24]. Hence in this case the bound on rrih can 
be satisfied even for relatively small mo and small mi/2- This revives the region where LSP 
pair annihilation is significant, which would otherwise remain forbidden for Aq = { see, 
e.g., LEPSUSY working group figures in Ref.[25]) ^. Moreover the low mo — mi/2 regions 
of the mSUGRA parameter space are characterized by relatively light squarks and gluino. 
Thus this region will be extensively probed at the early stages of the LHC experiments. It is 
therefore worthwhile to analyze the anticipated collider signals corresponding to this region 
in detail. 

Even for moderate negative values of Aq the WMAP allowed parameter space extends 
considerably. More importantly a small but interesting region where LSP pair annihilation 

^We follow the standard sign convention of Rcf. [23] for the signs of /i and ^o- 

^For Aq > 0, one requires mo and mi^2 typically larger than the corresponding values for Ao=0. This 
does not lead to any novel collider signal. 
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produces an acceptable dark matter relic density is revived. There is also a region where LSP 
- fi coannihilation is still the most important mechanism for having observed dark matter 
of the universe. Remarkably, even this region corresponds to a much smaller mi/2 compared 
to what one would obtain for the Aq =0 case. Consequently the squarks and gluinos become 
relatively lighter. Furthermore large negative values of Aq leads to a relatively light top 
squark with obvious characteristics in the collider signals. In fact for a small region of the 
parameter space, the LSP- ti coannihilation [12] may significantly contribute to the observed 
dark matter density. 

In this paper emphasis will be given on the features of the sparticle spectrum and signals 
at the Tevatron and the LHC corresponding to the WMAP allowed regions of the parameter 
space opened up by non-zero Aq. Such signals will be compared and contrasted with the 
expectations from the scenarios with Aq = 0. 

Several earlier dark matter analyses considered non-zero trilinear couplings in various 
SUSY models including mSUGRA[26]. More recently it has also been noted in the literature 
[27] that the DM allowed mSUGRA parameter space is expanded for non-zero Aq inspite of 
the constraint in Eq.l and the lower bound on rrih from LEP2. However, in this work we 
go further and point out the dominant annihilation/co-annihilation mechanisms that would 
produce the acceptable amount of neutralino relic density in different regions and analyze 
the novel collider signals associated with them. 

We have mainly considered the direct constraints from LEP2 searches and WMAP data 
on the mSUGRA parameter space. Other theoretical and experimental indirect constraints 
have also been considered in the literature. For example, it is well known that large values 
of Aq may lead to a charge and colour breaking (CCB) minima[28]. We shall comment on it 
in the next section. 

We note in passing that the constraints from the observed branching ratio BR(6— > 
57) [29] disfavours the mSUGRA parameter space where both mo and mi/2 happen to be 
small[22]. Parts of parameter space at the focus of attention of this paper belongs to this 
category. However, the theoretical prediction of this branching ratio in mSUGRA is based on 
the assumption of perfect alignment of the squark and quark mass matrices. This essentially 
means that the mixing angle factors at certain vertices ( e.g., the quark-squark- gaugino ver- 
tices) are identical to the Cabibbo-Kobayashi-Maskawa (CKM) factors at the corresponding 
SM vartices. This assumption is unnecessarily restrictive and may be falsified by small off 
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diagonal terms in the squarks mass matrices at Mq which may change the mixing pattern 
in the squark sector drastically. It should be emphasized that the above small off diagonal 
entries do not affect the ffavour conserving processes like neutralino annihilation or events at 
hadron colliders in any significant way. For a brief review of the model dependent assump- 
tions in the h— > 57 analyses we refer the reader to Okumura et al in [30] and Djouadi et al 
in [22]. 

2 The DM Allowed Parameter Space for Non-zero Aq 
and the Sparticle Spectrum 

In this section we relax the rather arbitrary choice Aq = and reexamine the scenarios with 
small mo and small mi/2 consistent with the WMAP data. We remind ourselves that this 
particular region is strongly disfavored by the rrih bound from LEP2 for vanishing Aq. 

We list the non-trivial consequences of considering a no n- vanishing Aq. First, the lightest 
Higgs boson becomes heavier even if we choose a moderately negative value for Aq keeping all 
other mSUGRA input parameters fixed. This happens via the mixing effects in the radiative 
corrections to rrih [24]. Moderately large negative values like Aq ^ —500 GeV are adequate 
to satisfy the experimental lower bound on rrih even for small tuq and small mi/2- 

However, three important additional effects should be taken into account while exploring 
the consequences of non- vanishing Aq{< 0). i) With a negative Aq having moderately large 
magnitude (thereby At - the trilinear coupling at the top quark sector at the weak scale 
being further negative), the lighter top squark ti will become much lighter than the other 
squarks because of the off-diagonal element mt{At — yucot /3) in the top squark mass matrix. 
In fact ti may even become the next lightest super particle (NLSP) or even violate the 
current experimental lower limit if At is sufficiently large and negative. We note here that 
unless tan (3 is very small, the amount of mixing in the ti sector is dominated by At in spite 
of the fact that \fi\ determined by REWSB also becomes larger as Aq is driven towards large 
negative values, ii) In the f sector fi becomes lighter for large \Aq\ for entirely different 
reason. Here /itan/3 dominates over Ar even for large |y4o| and moderate tan /3. This 
pushes fi to become lighter. The effect of mixing is obviously larger for larger tan /5 or a 
larger mi/2- As a result fi and to a lesser extent ti may become NLSP in different zones 
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of the parameter space around the bulk annihilation region. These two sparticles may thus 
coannihilate with the LSP and consistency with the WMAP limits on relic density may be 
obtained, iii) As mentioned in the introduction when |y4o| becomes large one may hit a CCB 
minima of the scalar potential [28]. We have imposed the CCB conditions at the electroweak 
scale in this analysis. However these constraints can be altogether relaxed if it is assumed 
that the universe is built on the EWSB false vacuum with a life time larger than the age of 
the universe [31]. We therefore do not pursue these constraints at other scales. 

In this analysis we have generated the electroweak scale sparticle spectrum from the 
input parameters at the unification scale by using the SuSpect code [32] with including the 
CCB condition as mentioned before. This code employs the two-loop renormalization group 
equations and implements the radiative electroweak symmetry breaking mechanism. We 
have used the code micrOMEGAs[33] for computing the neutralino relic density. We have 
included LEP2 lower limits for sparticle masses [20] and set top mass rrit =173 GeV. In 
particular the lower bound 94 GeV on m~ is quite important. However, we should also 
note here that the Higgs mass bound mh > 114 GeV should be treated with caution. The 
theoretical prediction for the Higgs mass in the MSSM has been computed including two loop 
corrections [24]. It is well known that this prediction involves an uncertainty of about 3 GeV 
[34, 35, 36, 37, 38]. The sources of this uncertainty are the momentum-independent two-loop 
corrections, the momentum-dependent two-loop corrections, the higher loop corrections from 
t/t sector etc. While presenting the allowed parameter space we have , therefore, delineated 
the regions corresponding to 111 GeV < rrih < 114 GeV. This parameter space which cannot 
be presently ruled out with certainty will be referred to as the uncertain rrih zone in this 
paper. 

Additionally, we have randomly varied rriQ by 5 GeV around the quoted value of 120, 80 
and 170 GeV in Fig.l, Fig.2 and Fig. 3 respectively. We have introduced this variation since 
different spectrum generators do have small uncertainties in computing the sparticle masses 
in the mSUGRA framework. Similarly we have used a variation of 10 GeV for in Fig. 4 
around Aq = -700 GeV. 

We have also shown the consequences of the indirect constraint from the measured 
branching ratio of 6 — * 57 [29, 30]. However as mentioned in the section. 1 these constraint 
becomes weaker if some additional theoretical assumptions are dispensed with. We have 
used the constraint 
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2.77 X 10"^ < Br{b 57) < 4.33 x 10"^ (3) 

at three sigma level. 

Fig. (1(a)) shows different WMAP allowed zones in the {Aq — mi/2) plane for tan/5 = 10, 
mo = 120 GeV. As mentioned before we have identified regions corresponding to 111 GeV < 
nih < 114 GeV (the uncertain zone) and rrih > 114 GeV (the regular zone). We find WMAP 
allowed regions for small mi/2 in both the zones. Additionally contours for m^-^ are also 
shown. The WMAP allowed parameter space has three distinct regions. 

1. In the region marked by the red dots the LSP pair annihilation is the dominant mech- 
anism. However, this alone cannot satisfy the relic density constraint. Additional 
contributions come from a significant degree of LSP-ti coannihilation. As already 
mentioned this coannihilation is a direct consequence of a top squark NLSP with mass 
as low as ~ 105 GeV due to large mixing effects in the top squark mass matrix for 
large and negative Aq. Such a light rrir with mr — m~o ~ 30 GeV is certainly within 
the striking range of Run II [39]. Thus if this scenario is indeed the one chosen by 
nature, the discovery of SUSY at LHC will be heralded by the discovery of the light 
top squark at the Tevatron. 

2. In the region marked by the blue (deep shaded) dots neutralino pair annihilation is the 
main mechanism for satisfying the observed relic density. However there are points, 
the ones with relatively large mi/2 in particular, where ~ 50% of the relic density is 
due to fi-LSP coannihilation. 

3. There is a third WMAP allowed region for larger mi/2(the pink/light shaded dots). 

Here fi-LSP coannihilation dominantes over LSP annihilation in satisfying the relic 
density. This coannihilation process becomes more and more effective as mi/2 in- 
creases. However as noted before, in a scenario with REWSB, a moderately large 
negative Aq yields large \fi\. This causes m^^ to become smaller for smaller values of 
mi/2 in comparison to the vanishing Aq case. Thus larger \Aq\ may trigger fi-LSP 
coannihilation at smaller values of mi/2 which in turn predicts lighter squarks and 
gluinos in the WMAP allowed region. 

The situation changes significantly if we change mo in either direction. Fig. (1(b)) illus- 
trates the case of mo = 80 GeV. Here |y4o| could not be as large as in the previous case 
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Figure 1: WMAP allowed region in the mi/2 — Aq plane for tan/5 = 10 and mo = 120, 80 
and 170 GeV. Lighter stop masses are shown in (blue) solid lines. Dashed lines correspond 
to lighter Higgs masses. Dot-dashed lines refer to 6 ^ 57 limits. WMAP allowed zones: 
The brown(medium shaded) region near zone (a) refers to LSP pair annihilation. In the red 
region near zone (b) LSP pair annihilation and LSP-ti coannihilation jointly produce the relic 
density. The blue(deep shaded) region near zone (c) refers to a mixture of bulk annihilation 
and LSP-fi coannihilation. The pink(light shaded) region near zone (d) corresponds to a 
dominantly LSP-fi coannihilation. Disallowed zones: Region (a) is disallowed because either 
xf falls below the LEP limits or rrih becomes unacceptably small (< 108 GeV). Region (b) 
is ruled out by the lower limits on mj^ or/and rrif^. Region (c) is disfavoured because CCB 
violating minima occurs at the elctroweak scale. Region (d) is discarded as fi becomes LSP. 
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Fig. (1(a)) because that would reduce m^^ below the experimental limit. This on the other 
hand does not allow any appreciable increase oirrih due to Aq. Smaller \Ao\ does not reduce 
m^^ appreciably so as to have any LSP-ti coannihilation. We note in passing that nif^ < 200 
GeV is disallowed since fi takes the charge of LSP here. This of course disfavours any LSP 
- ti coannihilation and ti is unlikely to be visible at the Tevatron [39]. 

The pure bulk region (shown in brown) is extended in comparison to the case of mo = 
120 GeV of Fig. (1(a)). However most of the extended parameter space corresponding to the 
bulk region has rrih < 111 GeV in this case. Additionally there is a sizable mixed region 
(where both LSP pair annihilation and LSP-fi coannihilation are important) which falls in 
the uncertain rrih zone. 

On increasing mo to 170 GeV (Fig. (1(c))) we find no bulk annihilation region. Here 
negative and larger Aq tends to reduce m^^, but a larger mo compensates. Thus ti may 
become the NLSP without having fi below Xi- As a consequence of such large Aq one finds 
region with rrih >114. This scenario has a large LSP-fi coannihilation region (pink/light 
shaded dots) as well as a small LSP-ti coannihilation region (red dots). However, part of 
the LSP-fi coannihilation region is associated with significantly light squark/gluino masses 
compared to the Aq = case and such regions are well within the reach of the early probes 
at LHC. 

Now we mention the status of the b— > s 7 constraint. The constraint disfavours the LSP- 
ti coannihilation region altogether in all the cases discussed above. But, segments of all three 
WMAP allowed regions are consistent with this constraint for mo = 80 GeV whereas a tiny 
part of mixed region and broad part of coannihilation region for mo = 120 GeV survive this 
constraint. On the other hand, for mo = 170 GeV the entire LSP- fi coannihilation region 
is allowed by this constraint. However, as mentioned in the section. 1 the 6 — 57 constraint 
may loose much of its impact if additional theoretical assumptions like the alignment of the 
quark and the squark mass matrices are relaxed. We therefore prefer not to exclude any 
region allowed by the WMAP data on the basis of this constraint. 

Fig. (2) shows the result for tan/3 = 5. Here, allowed values of Aq in Fig. (2) has a 
range larger than that of Fig.(l). The WMAP allowed regions for mo = 120 GeV Fig. (2 (a)) 
consist of only coannihilations of LSP with ti and fi in contrast to the Fig. (1(a)) where 
a significant amount of mixed zone also exist. The dark matter satisfied parameter space 
for mo = 80 GeV of Fig. (2(b)) survives only if full allowance is made for the theoretical 
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uncertainty in the lighter Higgs mass calculation i.e the WMAP allowed region coincides 
with the uncertain rrih region. However, most of the region satisfy the b ^ s'-f constraint of 
Eq.(3) except the pure bulk region and the LSP - ti coannihilation region. If this constraint 
is relaxed light top squarks within the striking range of Run II experiments are permitted. 
However, search for such light top squarks requires caution. For low tan/3 the branching 
ratio (BR) of the 4-body decay of this squark may be sizable or may even overwhelm the 
conventional loop decay [40]. Search strategies should be modified to include this usually 
neglected decay mode[41]. 
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Figure 2: Same as Fig.l except for tan/5 = 5. 



As an example of a large tan (3 case Fig. (3) shows the result for tan /? = 30 and mo = 
170 GeV. Here the WMAP constraint is satisfied via LSP-fi coannihilation even if nih > 114 
GeV. One also obtains a mixed region where pair annihilation is the dominant mechanism. 
However this region corresponds to the uncertain ruh region i,e 111 GeV < rrih < 114 GeV. 
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Figure 3: Same as Fig.l except for tan/3 = 30. 



We next study the available parameter space in the ttiq - mi/2 plane for a fixed negative 
value of Aq = —700 GeV. As we see from Fig. (4) we get a pure bulk annihilation region, a 
mixed region as well as a fi-LSP coannihilation region all satisfying the lighter Higgs bound 
and the WMAP limits. Additionally there is a significant region with LSP-ti coannihilation. 
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Figure 4: Same as Fig.l except for being drawn in mo — rni/2 plane for tan/3 = 10 and 
Ao = -700 GeV. 
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3 The Novel Collider Signals for non-zero values of the 
trilinear coupling 

In this section we simulate all possible squark - gluino events using Pythia (version 6.409) 
[42] at the LHC energy (y^=14 TeV). We compute the relevant BRs using the program 
SDECAY [32]. From Fig (1(a)) we have chosen the following benchmark points as shown in 
Table 1. 

A) This point corresponds to a mixed region where neutralino pair annihilation plays the 
dominant role. B)This point corresponds to the stau coannihilation region for a large negative 
value of Aq and C) This point corresponds to the lowest value of mi/2 in the LSP- fi 
coannihilation region for Aq = 0. 



mSUGRA 
parameters 


A 


B 


C 


mo 


120.0 


120.0 


120.0 


mi/2 


300.0 


350.0 


500.0 


Ao 


-930.0 


-930.0 


0.0 


tan 13 


10.0 


10.0 


10.0 


sgn{^) 


1.0 


1.0 


1.0 



Table 1: The three benchmark mSUGRA scenarios: (A): This parameter point with non- 
vanishing Aq belongs to a mixed region characterized by LSP pair annihilation and LSP- 
fi coannihilation, (B): This point with non- vanishing Aq belongs to a region where stau 
coannihilation dominates among the annihilation/coannihilation processes, (C) This point 
with vanishing Aq refers to the lowest value of mi/2 in the LSP-fi coannihilation region. All 
parameters with dimensions of mass are in GeV. 

The sparticle spectra for the three benchmark points are shown in Table 2. Table 3 
shows the branching ratios of the dominant decay modes of gluino and for the squarks of 
the first two generations. The branching ratios of the third generation of squarks are shown 
in Table 4. We have further computed the branching ratios of the dominant decay modes of 
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the lighter chargino and the second hghtest neutrahno. The result is shown in Table 5. 



S quark/ 
Slepton/Gluino/ 
Gaugino masses 


A 


B 


C 


9 


719.0 


827.0 


1149.0 


Ql 


670.0 


769.0 


1055.0 


Qb. 


647.0 


739.0 


1015.0 


h 


286.0 


396.0 


807.0 


h 


644.0 


728.0 


1016.0 


hi 


558.0 


652.0 


974.0 




639.0 


731.0 


1009.0 


k 


239.0 


268.0 


355.0 




226.0 


255.0 


346.0 


Ir 


169.0 


182.0 


224.0 




132.0 


148.0 


217.0 




218.0 


248.0 


345.0 


r2 


241.0 


268.0 


357.0 


xf 


232.0 


272.0 


386.0 


xt 


616.0 


669.0 


649.0 




121.0 


142.0 


205.0 


x'2 


232.0 


272.0 


386.0 


h 


117.0 


117.0 


115.0 



Table 2: The sparticle spectra in the three scenarios. All masses are in GeV. 

Next we present the total squark-gluino production cross-sections for the three scenarios 
in Table 6. These lowest order cross sections have been computed by CalcHEP (version 
2.3.7) [43]. From this table it is obvious that even for 10 fb~^ of integrated luminosity scenario 
A) will produce a remarkably large number of squark-gluino events compared to the other 
two scenarios. Thus it can be easily tested by the early LHC runs. On the otherhand if Xi 
- fi coannihilation happens to be the dominant mechanism for obtaining the present day 
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Decay modes 


A 


B 


C 


(squark/gluino) 


% 


% 


% 


9 qiq 


8.8 


10.4 


18.0 


9 QrQ 


18.4 


21.2 


35.2 


9 bib 


19.4 


18.6 


13.8 


9 hb 


5.6 


6.4 


9.4 


9 tit 


47.0 


43.0 


23.0 


Ql XiQ 


66.0 


66.0 


64.6 


Ql X2(l 


32.6 


32.8 


32.0 


Qr XiQ 


100.0 


100.0 


100.0 



Table 3: The BRs (%) of the dominant decay modes of the gluinos and the squarks belonging 
to the first two generations. 

thermal abundance and Aq is large negative (scenario B), the size of the signal may still be 
more than an order of magnitude larger than what is expected in scenario C) corresponding 
to = 0. The relative sizes of the total cross sections in the three cases can be qualitatively 
understood from the sparticle spectra as shown in Table 2. 

We next turn to some distinctive features of the signals from the three scenarios. This 
criteria may be utilized as a confirmation of the underlying SUSY scenarios that correspond 
to different mechanisms of neutralino annihilation for an acceptable relic density. From 
the decay tables it is evident that the signals from decaying sparticles in A) and B) will 
contain many more tau-leptons than electrons or muons showing thereby a strong departure 
from "lepton universality". However, in scenario C) this universality in the signal will hold 
approximately. 

The main reason for this "non-universality" lies in the dominant 2 - body decay modes 
of xf and X2- is to be noted that in scenario A) the xf decays into the lighter stau(with a 
dominant tr component) with a large BR. This is due to the fact that the lighter charginos 
can decay only into right handed sleptons via two body modes. However this decay is allowed 
only through the small higgsino component of the gaugino like chargino. Thus the decay 
into Ti dominates. Again since z/^ is lighter than the other sneutrinos, xf decays into Uj- 
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Decay modes 
(heavy squark) 


A 

% 


B 

% 


C 

% 


bi ^ x'b 


12.0 


14.0 


19.7 


bi Xit 


16.0 


21.5 


35.0 


bi tiW- 


70.0 


63.4 


13.0 


bi X2t 


0.0 


0.0 


29.0 


ii Xib 


100.0 


49.0 


36.0 


ti ^ x\t 


0.0 


51.0 


29.0 


ti ^ xlt 


0.0 


0.0 


14.0 


b2 ^ X'lb 


28.4 


3.3 


24.5 


b2 ^ Xlb 


3.9 


4.2 


6.7 


h Xit 


5.9 


6.9 


12.0 


1-2 ^ tiW~ 


61.0 


54.0 


13.0 



Table 4: The BRs (%) of the dominant decay modes of the squarks belonging to the third 
generation. 

and tau with a sizable BR. In all three scenarios the sneutrino decays invisibly via the mode 
i>j — > z/j + (i=e,/i, r). Thus the second lightest neutralino dominantly decays either into 
invisible neutrino-sneutrino pairs or into r and ti with large BRs. As a result there is a 
significantly large tau excess in the semi- inclusive signal r^+X,-, (where stands all final 
states without a tau) compared to e^+Xg or /i^'^+X^. 

In the squark gluino events xt ^"^^ X2 primarily arise from the decays of and ti 
decays into Xa ^iid b with 100% BR. Moreover since ti is much lighter compared to the 
other squarks due to the large magnitude of the Aq parameter the gluinos decay into ti — t 
pairs with a large BR. 

Scenario B) has all the above features albeit to a lesser extent. It is particularly important 
to note that more than 50% of the tis decays into t — x\ pairs. The decays of the t tend to 
restore lepton universality. Thus the excess of events with r leptons over the ones with e or 
/i is reduced to some extent. 

In scenario C) the excess of tau leptons is reduced drastically as is evident from decay 
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Decay modes 
(Gauginos) 


A 

% 


B 

% 


C 

% 




78.6 


46.0 


8.4 




5.2 


24.0 


33.8 




15.0 


24.0 


17.8 




2.6 


2.6 


5.1 




5.4 


25.0 


36.0 




16.0 


25.0 


19.2 


xt T^tVr 


76.0 


44.0 


7.9 



Table 5: The BRs (%) of the dominant decay modes of the hghter chargino and the second 
hghtest neutrahno. All sneutrinos decay into the invisible channel + Xi in the three cases 
under study. 

tables. 

However, it must be borne in mind that the decays of the standard model particles 
(mainly t, Z and W) present in the signal exhibit lepton universality. Moreover, the heavy 
flavour (6, c, t) decays yield more electrons and muons than taus. It is therefore essential to 
calculate the relative abundance of different charged leptons in the signal through simplified 
simulations. 

To begin with we ignore the difference in the detection efficiencies of different leptons. 
This will be addressed later in the paper. Again the production and decay of all squark- 
gluino pairs are generated by Pythia. In this simplified parton level analysis all SM particles 
other than the Z and t are treated as stable. All unstable sparticles are allowed to decay. 
Hadronization, fragmentation and jet formation are switched off. We also impose the cuts 
Pt > 30.0 GeV and {rjil < 2.5 on e,/i and r. This is to have a feeling for the numbers when 
e's and the /i's in the final state from b and c decays are minimized (see the improved analysis 
presented below). The results are shown in Table 7. In this table the first row corresponds 
to the number of semi-inclusive events with only one r. The second and third rows give the 
corresponding numbers for final states with one e and one /i respectively. 

From Table 7 it is obvious that in scenario A) the number of events with only one r is 
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Process 


a{ pb) 


A 


B 


C 




1 41 


55 


05 

\j .\JO 




2.35 


1.08 


0.15 






1 1 8 


n 1 (S 




u. / o 




n 1 




n 18 


n 08 


n ni 


HRhR 






n 1 1 






1 


01 


(It fin 
HLHR 


0.47 


0.25 


0.05 


Qt On + C.C 


0.43 


0.21 


0.04 


tin 


6.12 


1.14 


0.02 




0.07 


0.03 


0.004 


tit^ + C.C 


0.001 


0.0002 


1.5x10"^ 




0.17 


0.07 


0.005 




0.08 


0.03 


0.004 


Total 


15.58 


5.74 


0.74 



Table 6: The production cross sections of all squark-gluino events studied in this paper. 

much larger than the corresponding events with e or /i. The excess of events involving only 
one T lepton is also seen in scenario B). The scenario C) exhibits lepton universality. 

Since the gluinos dominantly decay into b ot t flavoured squarks in scenarios A) and B) 
we expect a sizable number of b quarks in the final state. This is also illustrated in Table 7 
(see the last five rows). 

Thus in principle 6-tagging may be effectively used to suppress the backgrounds except 
the ones from t — i production. 

To demonstrate that the above tau excess survives the typical kinematical cuts designed 
for SUSY search experiments at the LHC we proceed as follows. Since a detailed background 
simulation is beyond the scope of this paper we apply a slightly modified version of the generic 
cuts introduced by the ATLAS collaboration to eliminate the SM backgrounds in their study 
of inclusive SUSY signals [44]. We shall, however, study the response of the t — t events. 
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A 


B 


C 


iT + Xr 


29870 


11860 


1340 




5274 


4251 


1260 




5294 


4232 


1262 


Ir + 06 + 


10750 


4581 


747 


Ir + 16 + 


19 


13 


1 


Ir + 26 + Xr 


17099 


6589 


507 


Ir + 36 + 


9 


6 





Ir + 46 + Xr 


425 


665 


80 



Table 7: The number of semi-inclusive events with one r, one e and one /i at the parton level. 
X stands for all possible final states excluding the lepton indicated by the subscript. The 
last five rows indicate the number of 6 quarks in the r type events. All numbers correspond 
to an integrated luminosity of 10 fb^^. 





A 


B 


C 


ti 


Ir + X, 


3113 


1402 


239 


481 




1179 


1246 


820 


1295 


le + Xe 


1138 


1263 


829 


1354 



Table 8: The number of semi-inclusive events with one detected r jet, one isolated e and one 
isolated /i computed by Pythia. The selection criteria and the kinematical cuts are given the 
text. The last column gives the contributions of the background from t — i production. 

likely to be the dominant SM background, to these cuts after our signal analysis. 

Again all squark-gluino events are generated by Pythia. Initial and final state radia- 
tion, decay, hadronization, fragmentation and jet formation are implemented following the 
standard procedures in Pythia. We impose the following selection and background rejection 
criterion: 

1. Only jets having transverse energy > 30 GeV, pseudo rapidity |?7j| < 4.5 and the 
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jet-jet isolation parameter Ai?(ji,j2) > 0.5 are selected from the toy calorimeter of 
Pythia. 

2. The jets arising from hadronic tau decays are selected as follows. Firstly r's having 
transverse momentum pt > 20 GeV and pseudo rapidity |?7t-| < 3 are identified. These 
rs are matched with the jets by setting the following criterion AR{T,j) < 0.4 and 
Ej^/E^ > 0.8, where Ej, is the transverse energy of jets and E^ is the same for the r 
being matched. The possibility of detecting the matched rs is then computed using 
the detection efficiencies in [45](see Fig 12.9, p 444 in section 12.1.2.2). 

3. We also require that an isolated e or /i in an event have px > 30 GeV and the lepton-jet 
isolation parameter AR{l,j) > 0.5, where / stands for either e or /i. The detection 
efficiencies of these leptons are assumed to be 100%. 

4. We reject events without at least two jets having transverse momentum > 150 GeV. 

5. Events are rejected if missing transverse energy < 200 GeV. 

6. Only events with jets having St > 0.2, where St is a standard function of the eigen 
values of the sphericity tensor, are accepted. 

The strong pt and isolation criteria on e and fi enable us to exclude the most of leptons 
from semileptonic b and c decays from the list of isolated leptons in an event. We have 
checked that a suitable cut on M^ff, where M^fj = |^| + Ilf^i\pj.\ + Tif^^\pl^\ {I = e, /i) does 
not improve the ratio S/ \/B significantly, where S corresponds to number of signal events 
and B represents the number of background events. 

We present the number of semi-inclusive events with only one detected r jet, with only 
one isolated /i or e in Table 8. The number of events corresponds to an integrated luminosity 
of 10 fb~^. It is clear from the table that inspite of non-ideal r detection efficiencies there is 
indeed an excess of events with tau leptons over that with electrons and muons in scenario 
A). In scenario B) the number of events with tau leptons is reduced as expected compared 
to A) but is still significantly larger than the corresponding number in C). 

We note in passing that there are several proposals for detecting the taus [45]. The 
combined efficiency may be larger than the efficiency obtained from tracker isolation alone 
which we have used. However, the combined efficiency is process dependent [45]and in this 
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simplified analysis we have not employed it. Also following the CMS guidelines we have 
not considered the possibility of detecting tau-jets with < 30 GeV. However, we have 
noted that a large number of tau-jets in the final state have pt in this region. If future works 
establish the techniques of detecting these relatively low pt tau jets, the excess of final states 
involving rs will be even more dramatic. 

We next study the dominant background from t — i events using the selection criteria 
and kinematical cuts given above. The leading order cross-section for this process computed 
by CalcHEP is 400 pb. 

The results are given in the last column of Table 8. The S/\fB ratios in the three sce- 
narios indicate that statistically significant signals are expected in all cases. Since 6-tagging 
is not particularly helpful in rejecting this particular background we have not incorporated 
it here. However in a more complete background analysis 6-tagging may indeed be a useful 
tool as our parton level analysis indicates. 

As noted in the introduction a consequence of large negative Aq is the possibility of 
a light top squark which can be as light as the 105 GeV (Fig.l - 4). For such a light ti 
its coannihilation with the LSP may be one of the important mechanisms for having the 
observed DM in the universe. Moreover, the LSP mass in this region of the parameter space 
is about 25-30 GeV lesser than m^^. Hence, the top squark NLSP will decay via the loop 
mediated mode c+x\ with 100 % BR. The prospect of light stop squark search in this channel 
at Tevatron Run II was studied in [39] . If this scenario for thermal abundance of dark matter 
in the universe is realized in nature SUSY search at LHC will certainly be heralded by the 
discovery of ti at the Tevatron. Even in scenario B) the ti will be considerably lighter 
compared to scenario C) (see Table 2) although its mass will still be outside the kinematic 
reach of the Tevatron. However this squark will have a sizable production cross section at the 
LHC. Hence its discovery may be an additional confirmation of a mechanism of producing 
an acceptable DM relic density based on x5 - coannihilation with large negative Aq. 

An additional feature of WMAP allowed regions of the parameter space with low mo is 
the existence of sneutrinos decaying invisibly. This happens in all three scenarios A, B and 
C. Moreover, the x\ decays into sneutrino neutrino pairs with sizable BR. Thus in addition 
to the LSP there are two other carriers of missing energy. An important consequence of these 
scenarios is the spectacular enhancement of the signal e^e~ — > 7 + missing energy over the 
SM background from e^e~ — > ViVi + 7 (i=e,/i,r). The signal comes from three channels 
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e^e~^ XiXi7;X2X275 The signal from the first SUSY process was computed in [46] 

with the approximation that Xi is a pure photino. All three cross-sections including the 
effect of neutralino mixing were computed in [47]. The results of [47] were compared with 
that of CalcHEP in [48]. The observation of the signal at the ILC may provide additional 
confirmation of SUSY dark matter if produced in the low mo, mi/2 region of the mSUGRA 
model. 

4 Conclusions 

In this work we observe that parts of the low itiq - mi/2 region of the mSUGRA model with 
moderate or large negative values of the common trilinear coupling Aq at Mq are consistent 
with the current direct bounds on sparticle masses from LEP as well as the WMAP bound 
on the dark matter relic density (Fig.l - 4). In particular the LEP bound rrih >114 GeV 
can be easily satisfied for relatively small mo and mi/2 if is large and negative (Fig.l - 
4). This possibility is excluded by the adhoc choice Aq =0 made in many current analyses. 

The relic density is produced by several processes including LSP - pair annihilation, LSP 
- ti or LSP - fi coannihilation or by suitable combinations of these processes. 

This parameter space corresponding to relatively light squarks and light gluinos is im- 
portant since it will be extensively scanned at the early runs of the LHC. A large excess of 
semi-inclusive events of the type r'^+Xr (see Tables 7 and 8) over e^+Xg or /i^+X^ events 
is a hallmark of scenarios with large negative Aq in mSUGRA. Our simulations with three 
bench-mark scenarios (see Table 1) consistent with direct lower bounds on sparticle masses 
and WMAP data establish this for a modest integrated luminosity of 10 fb~^. 

A natural consequence of large negative Aq is a light ti. If that be the case then ti - Xi 
coannihilation may be an important mechanism for the present day dark matter density. In 
this scenario the discovery of SUSY at the LHC is likely to be preceded by the discovery of 
ti at the Tevatron. 

Several groups have also imposed indirect constraints on the allowed parameter space 
[22]. However, it is well known that each of these indirect constraints employ additional 
theoretical assumptions which are not fool proof. For example the requirement that there 
be no CCB minima of the scalar potential deeper than the EWSB vacuum [28] becomes 
redundant if the latter is assumed to be a false vacuum with a life time larger than the 
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age of the universe [31]. Similarly the constraint from the measured branching ratio of the 
decay b ^ s'-f derived under the purely theoretical assumption that quark and squark mass 
matrices are perfectly aligned can also be evaded if this assumption is relaxed ( see,e.g., 
Djouadi et al in [22]and references therein). Although we have shown the consequences of 
some of these indirect bounds in our figures it may not be prudent to discard any parameter 
space if these bounds are violated. 
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